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Device and method for forming an image composed of a plurality of sub-areas 



The invention relates to a device for forming an image which is composed of a 
plurality of sub- areas, with a detector which includes a plurality of sensor elements for 
generating image data and read-out units which are associated with the sub-areas of the 
image. The invention relates to a method of forming an image which is composed of a 
5 plurality of sub-areas, with each sub-area there being associated a respective read-out unit 
and image data of adjoining image areas of neighboring sub-areas being evaluated in order to 
y mitigate differences in amplification characteristics. The invention also relates to a radiation 

Si examination device for performing the method according to the invention which includes a 

radiation source and a detector device with an appropriate device for carrying out the method. 
P 10 The invention also relates to a computerprogranmi for the correction of image data of an 
U image which is composed of a plurality of sub-areas. 

Methods and devices of this kind are used, for example in image generating 
devices of X-ray or other radiation examination devices. 

WO 96/19893 discloses an image pick-up apparatus in which an image signal 
15 is decomposed into a plurality of sub-images. More specifically, the image signal is 

decomposed into two sub-images which are applied to respective detector fields. An overall 
image signal is formed from these sub-images again in a combination unit. A correction unit 
is provided so as to correct differences in the sub-images. To this end, brightness values of 
columns of a first sub-image are multiplied by a column amplification factor and brightness 
20 values of rows of a second sub-image are multiplied by a row amplification factor. The 
amplification factors are derived from reference images and are stored in a memory. 

For most fields of application of said methods and devices, notably for 
medical X-ray detectors but also for ultrasound apparatus, computed tomography apparatus 
or magnetic resonance apparatus, it is very important that the images formed are artefact-free 
25 to a high degree. Medical images are commonly reproduced as grey scale images. For 

example. X-ray and ultrasound apparatus or computed tomography or magnetic resonance 
tomography apparatus are used to form such grey scale images. 

The invention proposes a special device for a flat dynamic X-ray detector for 
the correction of images that contain defects. Such X-ray detectors are used for X-ray 
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examination apparatus in the medical diagnostic field. They can be used as universal detector 
components in various application-specific X-ray apparatus. For their operation in the various 
types of application it is very important that the images produced are free from artefacts to a 
high degree. 

5 Dynamic X-ray detectors are subdivided into a plurality of sub-areas, a 

respective read-out unit being used for each sub-area. Each sub-area comprises a plurality of 
image areas. The detected signals or image data of a sub-area are amplified by a read-out 
unit; one read-out unit then reads the plurality of image areas of a sub-area. Neighboring 
image areas that are situated at the edge of adjoining sub-areas, therefore, are read out by 
10 different read-out units. The inherently different amplification behavior of such different 
read-out units gives rise to abrupt grey scale transitions between such neighboring image 
S areas. Such abrupt grey scale transitions appear as stripes in the image to be formed. 

C| Such influencing is further intensified by noise suppression or image 

!t subtraction. The described artefacts become less visible as the noise component is greater. 

y I 

15 Because the human eye has an integrating property, the effectively present noise component 
12 is suppressed and the artefacts or defects become clearly visible. In the case of subtraction, 

required for the highlighting of image contents, two images are subtracted from one another, 
pj Such subtraction makes a difference in the amplification behavior of adjoining read-out units 

[T particularly evident. 

P 20 Such differences in non-linear amplification behavior can be corrected by 

means of regular recalibration processes and the use of the information concerning the non- 
linear amplification behavior of neighboring read-out units, which is derived during the 
recalibration processes. However, such a recalibration necessitates interventions by and - 
instruction of the medical or technical staff. 

25 In order to enable recalibration to be carried out, a predetermined number of 

X-ray exposures is performed with reproducible radiation doses incident on the X-ray 
detector. The amount of the radiation dose each detector sub-area can also be determined by 
the reproduction of special calibration phantoms, so that a deviation between the non-linear 
amplification characteristics can be derived from the differences occurring in the image. 

30 Because such non-linear amplification behavior is subject to change, 

correction must be performed regularly. The change cannot be considered as a fixed 
parameter over the service life of an imaging system. Such a non-linear amplification 
behavior can be induced by changes of the detector position, by inhomogeneous heat losses 
or also be the result of aging of the amplifier circuits. 




•PHD 99.198 

3 28.09.2000 
Such corrections that have to be performed manually at regular intervals have 
the drawback that on the one hand the attendant staff must intervene and that on the other 
hand additional X-ray exposures are required for the regular manual corrections. An 
adaptation table with correction data based on such manual correction would be used for all 
5 subsequent image data until a manual correction is performed again. Consequently, changes 
in the amplification behavior which occur in the mean time are not taken into account. 

Therefore, it is an object of the invention to provide a device and a method 
that enable regular and accurate recalibration or correction of non-linear amplification 
behavior to be performed without additional X-ray doses or user interventions being required. 
10 This object is achieved by means of a device for forming an image which is 

composed of a plurality of sub-areas, which device includes 
Q a detector which includes a plurality of sensor elements for generating image data, 

sj read-out units which are associated with the sub-areas of the image, 

an analysis unit which is arranged to evaluate image data from adjoining image areas of 
01 15 neighboring sub-areas and to generate correction data, and 

ru 

M= a correction unit which is arranged to correct incorrect image data by means of correction 

data. 

ry The X-ray detector is composed of a matrix of sensor elements which are 

lI arranged in rows and columns. The image produced by the X-ray detector can be subdivided 

y 20 into sub-areas, each sub-area including a plurality of image areas and a respective read-out 

unit being associated with each sub-area. The read-out units serve to amplify detected image 

data. 

The read-out units are constructed as amplifiers and each amplifier amplifies 
the image data applied thereto with an inherent non-linear amplification characteristic; this 
25 leads to deviating amplification of the image data. 

The invention is based on the idea that the image data of adjoining image areas 
of neighboring read-out units do not exhibit abrupt changes of value. It may indeed be 
assumed that the image data of neighboring image areas assume approximately the same 
values. Conversely, it can be deduced that the occurrence of an abrupt change of value may 
30 be attributed to non-linearities. 

The image data of medical images concern mainly grey scale values. An 
image area may be formed by an image column or image row or parts thereof and contains a 
plurality of pixels. It is assumed that the grey scale values of neighboring image areas have 
the same grey value scales. 
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For color images, in which the image data is represented by a vector 
containing the RGB components, it may also be assumed that adjoining image areas of 
neighboring sub-areas normally assume the same values. The same holds for images 
represented by amplitude values and phase values. 
5 In order to achieve high-quality correction, a series of special rules is formed 

on the basis of the described assumption. 

Pairs of grey scale values which belong to adjoining image areas of 
neighboring sub-areas and whose difference in grey scale value exceeds a selectable 
threshold value are not taken into account for the formation of the correction data, because it 
10 is very likely that this difference is not due to the non-linear amplification behavior of the 
neighboring read-out units. Furthermore, image areas in which an interpolation was 
performed by a preprocessing operation on the basis of pixel errors are also excluded from 
the formation of the correction data. Image areas which are situated in the direct vicinity of 
collimator shutters whose position is calculated by means of an algorithm are also not taken 
15 into account for the formation of the correction data. 

The grey scale values of neighboring image areas which are amplified by 
different amplifiers are analyzed in an analysis unit. The analysis unit produces correction 
data on the basis of image data received from adjoining image areas of neighboring read-out 
units. The correction data is applied to the incorrect image data in a correction unit in such a 
20 manner that the image data of, for example, the right-hand sub-area, amplified by means of 
the respective amplification characteristic, are adapted to the first amplification characteristic 
of the neighboring left-hand sub-area. This correction is then applied to all image areas 
belonging to this sub-area. 

The correction or adaptation is performed for all sub-areas, so that image 
25 defects due to different non-linear amplification characteristics no longer occur after this 
correction. 

The correction data is preferably stored in a memory. It can be fetched 
therefrom at any time so as to be used in the correction unit. 

Preferably, the image data generated by the detector is applied to the analysis 
30 unit only at a reduced clock rate. The calculation of the correction data requires an effort 

proportional to the required accuracy, so that, for example, only every second image from the 
detector is applied to the analysis unit. The accuracy of the correction is not degraded by this 
delay, because changes in the amplification behavior of the read-out units do not occur in 
these time relations but can be observed over prolonged periods of time only. 
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The analysis unit includes a histogram generator. The histogram generator 
receives image data of neighboring image areas of adjoining sub-areas. 

A respective histogram is formed concerning the occurrence of the 
corresponding grey scale values of each of the two image areas observed. Depending on the 
5 desired accuracy of the correction, the grey scale values are classified in classes. The relevant 
image data is applied to the histogram generator over a given period of time. These two 
histograms, in which the occurrence of the corresponding grey scale values of the individual 
image areas is reproduced over the entire range of feasible grey scale values, are applied to 
the summing unit. 

10 The sunrmiing unit forms a cumulative histogram each time from the 

histograms produced in the histogram generator. The resultant cumulative histograms also 

0 exhibit differences. A different amplification behavior of the read-out units amplifying the 
\i neighboring sub-areas can be attributed to these differences. 

^ From an analytical point of view the histogram generator forms density 

01 15 functions which are integrated in the summing unit, so that in this case the distribution 
u functions are formed from the density functions. 

The adaptation unit calculates a functional dependency of the image data of 
the second image area on the image data of the first sub-area. This functional dependency is 
taken up in an adaptation table as a confrontation of corresponding grey scale values. The 
20 correction data represent the contents of the adaptation table applied to the incorrect image 
data by the correction unit. 

The image data is corrected in the correction unit by application of an 
adaptation table of this kind in such a manner that defects or artefacts which are due to a , 
different amplification behavior are no longer visible in the X-ray image after this correction. 
25 A display unit displays a medical image, for example intended for diagnostic 

purposes, which has been freed from such errors. 

The calculation of the correction data in the analysis unit is not restricted to 
the application in the X-ray field, but can be used in all fields of signal processing in which a 
time-invariant non-linear amplifier behavior occurs. 
30 For the generation of the correction data it is also possible to use, in 

dependence on the application, a plurality of image areas of the sub-areas to be adapted, the 
density functions or histograms thereof then being linked by means of a suitable operation. 

In a further prefered embodiment, according to the invention the pixels 
situated at the boundary are observed and, while utilizing the image data of the adjoining 
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pixels of the neighboring sub-area, estimated values or reference values are determined for 
these pixels. These reference values indicate the value that the image value associated with 
the pixel should have if the read-out units were to have approximately the same amplifier 
characteristics. Thus, a correction value is determined by comparison of the actual image 
5 value of the relevant pixel and the estimated value determined; this correction value serves to 
correct the respective image value occurring in the sub-area to be corrected. 

In the simplest case the correction value may be a correction factor which is 
formed from the ratio of the estimated value to the image value. For the relevant image value 
and for the sub- area to be corrected this correction value can be taken up in an adaptation 
10 table or a so-called Look-Up Table. Needless to say that it is also possible to take up the 
product of the image value and the correction value, or in the above case the relevant 
estimated value itself, directly in such an adaptation table for the relevant image value. 

In the simplest case the image value of a directly adjacent pixel of the 
neighboring sub-area can be used as the estimated value. In this case it is assumed that no 
15 changes whatsoever occur at the boundary between one sub-area and the other. Such a 
correction, based merely on a comparison of individual pixels to each side of a boundary 
between sub-areas, does not require a large amount of calculation work and, therefore, can be 
quickly and simply carried out. However, it has the drawback that a slight gradient which 
actually occurs in the image and extends across the boundary is locally eliminated at the area 
20 of the boundary by calculation, and hence an error is produced again. 

In a further, preferred method which takes into account such gradients at the 
area of the boundary image values of pixels of the adjoining image area of the neighboring 
sub-area are extrapolated across the boundary in order to determine the estimated value. The 
term "adjoining image area", therefore, is not restricted to the directly adjoining row or 
25 column of the neighboring sub-area, for example, in the case of a detector consisting of rows 
and columns, but also covers a combination of several rows or columns or parts of rows or 
columns when this combination overall adjoins another sub-area by way of one side. 

In the simplest case this extrapolation is performed over only two pixels of a 
sub-area which are situated adjacently and transversely of the boundary line. Needless to say, 
30 however, that it is also possible to perform such an extrapolation over an arbitrary number of 
pixels which are adjacently situated pixels on one line. 

In an especially preferred embodiment a first correction value is formed first 
of all for the image value of a pixel of the sub-area to be corrected. Subsequently, for a 
neighboring pixel of the neighboring sub-area which directly adjoins this pixel of the sub- 
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area to be corrected there is determined an estimated value for the neighboring pixel while 
utilizing image values of the sub-area to be corrected. This is done in the same way as the 
determination of the estimated value for the pixel of the sub-area to be corrected itself. This 
means that an extrapolation over the values of the pixels of the neighboring sub-area is 
5 performed across the boundary for the pixel of the sub-area to be corrected; conversely, an 
extrapolation is performed over the image values of the pixels of the sub-area to be corrected 
in order to obtain the estimated value for the neighboring pixel. A second correction value is 
then formed by comparison of this estimated value for the neighboring pixel with the actual 
image value of the neighboring pixel. Finally, from these two correction values there is 
10 derived a common correction value which is applied to the relevant image values of the sub- 
area to be corrected. This common correction value may be the mean value of the two 
p correction values. A further possibility consists in that, for example, a weighted mean value 

is calculated. 

=p The estimation process can in principle be performed for all pixels situated 

g-l 15 along a boundary. In that case it is to be assumed that many of the pixels have the same 
rf image value. It would make more sense to form a common correction value for the relevant 

= image value in the sub-area to be corrected from the correction values of different pixels of 

hj the sub-area to be corrected which have the same image value. This can be performed, for 

\^ example by forming the mean value, the modulus or the median. 

Q 20 The correction values determined for the individual image values of the 

individual sub-areas are stored in an adaptation table and are fetched therefrom for the 
purpose of correction. The correction values may in principle again be correction factors 
which are applied to the relevant image value. However, they may, of course, also be 
previously corrected finished values simply substituted for the relevant image value. 
25 Because it is not likely that all possible image values occur along a boundary 

of a sub-area, for some of the image values this correction method does not yield data for all 
possible image values during the formation of a relevant adaptation table. However, any gaps 
arising can be readily filled in by interpolation between the correction values found by means 
of said method. In order to keep the required amount of calculation work small, in given 
30 circumstances it will even suffice to determine a correction factor directly for single image 
values only and to determine the remainder of the table by interpolation while spending only 
a small amount of arithmetic work. 

In conformity with the foregoing, the device for carrying out the method 
includes a memory for the storage of such an adaptation table. Using the adaptation table 
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stored in this memory, the image values are then always corrected in a correction unit for 
being displayed on a display unit. 

It is in principle possible to apply the correction values acquired by means of a 
given image directly to this image when the image data is delayed on its way to the 
5 correction unit. In order to increase the processing speed, however, it makes sense to use the 
correction values produced for previous images, or an adaptation table formed therefrom, and 
to check and correct this adaptation table in the memory each time by means of the new data. 

For example, a sliding mean value can then be generated over the contents of 
the adaptation tables acquired from successively formed images. This means that the 
10 correction values are averaged each time over the last images and it is useful to perform a 

weighting operation during the formation of the mean value. In that case it is advantageous to 

0 attach less weight to a new set of correction values than to the contents of the adaptation table 
sj present in the memory, because this table has already been formed from a mean value of the 

preceding images. Short-term fluctuations which could be due to other errors are thus 

y I 

01 15 suppressed. 

LI Because changes in the amplification behavior of the read-out unit usually do 

^ ^ not occur in time relations corresponding to the intervals between individual images of a 

rj series of images but can rather be observed over prolonged periods of time, the amount of 

12 calculation work required can in principle also be reduced by applying images to the analysis 

2 20 device only with given intervals in time, for example only every second or third image, in 

order to check or correct the adaptation table in the memory. 

The calculation of the correction data by means of said method is not restricted 
to applications in the X-ray field. It can be applied in all fields of signal processing wherein h 
a time-invariant, non-linear amplifier behavior occurs between different image areas of a 

25 common image. 

A radiation examination device according to the invention, notably an X-ray 
examination device, includes a radiation source and a detector device. The detector device 
includes a detector with a pluraUty of sensor elements, several sensor elements being 
combined so as to form sub-areas and each of said sub-areas being provided with a respective 

30 read-out unit. The detector device also includes an analysis unit for the evaluation of image 
data of adjoining image areas of the neighboring sub-areas and for producing correction data, 
the analysis unit including means for forming the estimated value by means of the method 
according to the invention and means for forming a correction value while utilizing the 
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estimated value. The detector device also includes a correction unit which serves to correct 
incorrect image data while utilizing the correction data acquired. 

The object of the invention is also achieved by means of a method of forming 
an image composed of a plurahty of sub-areas, wherein a read-out unit is associated with 
each sub-area and image data of adjoining image areas of neighboring sub-areas are 
evaluated in order to mitigate differences in amplifier characteristics. 

The object of the invention is also achieved by means of an X-ray examination 
apparatus which includes an X-ray source for emitting X-rays and for forming an X-ray 
image, an X-ray detector which serves to derive an optical image from the X-ray image, 
contains sensor elements which are arranged in rows and columns, and at least two associated 
amplifiers for reading out detected image data, at least one amplifier for reading out detected 
image data being associated with each sub-area, an analysis unit for generating correction 
data on the basis of the evaluation of image data of adjoining sub-areas of neighboring 
amplifiers, and a correction unit for the correction of the incorrect image data by means of 
the correction data. 

The object of the invention is also achieved by means of a computer program 
for the correction of image data of an image composed of a plurality of sub-areas, a 
respective read-out unit being associated with the sub-areas of the image and image data of 
adjoining sub-areas of neighboring read-out units being evaluated by the formation of a 
histogram in order to generate correction data after integration of the histograms, which 
correction data is used to adapt the image data of one sub-area to the amplifier characteristic 
of the read-out unit amplifying the adjoining sub-area. 



Further details and advantages of the invention will become apparent from the 
dependent claims and the following description in which an embodiment of the invention will 
be described in detail with reference to the Figures. Therein: 

Fig, 1 shows diagranmiatically the construction of a device for forming 

images, 

Fig. 2 shows the construction of a detector, 

Fig. 3 shows density functions of image data of adjoining image areas of 
neighboring read-out units. 

Fig. 4 shows probability functions of the image data. 
Fig. 5 illustrates the adaptation of the image data. 
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Fig. 6 shows an X-ray image containing defects, and 
Fig. 7 shows an X-ray image wherefrom defects have been removed. 
Fig. 8 illustrates diagrammatically the extrapolation for obtaining an estimated 
value at a boundary between two sub-areas; 
5 Fig. 9 shows diagrammatically the construction of a device for the formation 

of images; 

Fig. 10a shows an image of an X-ray detector subdivided into a plurality of 
sub-areas, with each pixel of a given sub-area there being associated an unambiguous grey 
scale value which is characteristic of the sub-area; 
10 Fig. 10b shows a distorted X-ray image of the detector shown in Fig. 10a; 

Fig. 10c shows the adaptation table in the form of grey scale values for the 
p correction of the image shown in Fig. 10b; 

J=[ Fig. lOd shows the distorted X-ray image of Fig. 4b after correction by means 

=P of an adaptation table as shown in Fig. 10c. 



15 



Fig. 1 shows the construction of a device according to the invention. The detector 1 1 
generates the image data which is applied to a preprocessing unit 19 and therefrom to the 
analysis unit 12 and the correction unit 13. For example, an offset correction is performed 
Q 20 already in the preprocessing unit. The histogram generator 15 in the analysis unit 12 
generates a respective histogram from the image data of adjoining image areas of 
neighboring sub-areas. The histograms are applied to the sunmiing unit 16 in which the 
cumulative histograms are formed from these histograms. The adaptation unit 17 calculates 
the corresponding grey scale value of the other image area for each grey scale value observed 
25 for the one image area. The functional dependency of the image data observed is produced on 
the basis of this calculation. The correction data thus produced is applied directly to the 
correction unit 13 and is also stored in a memory 14 in always the most updated form. The 
image data from the detector, applied directly to the correction unit 13, is then corrected by 
means of the correction data in the adaptation table and is transferred to the display unit 10 
30 which need not necessarily be attached directly to this device, but may also be present in a 
remote station; it is then connected thereto via a data network. 

Fig. 2 shows diagrammatically the typical construction of an X-ray detector. 
An X-ray detector of this kind converts incident X-rays into, for example, visible light which 
is detected by photosensors and applied to the read-out units. Also known are X-ray detectors 
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which convert the X-rays directly into electric charges. The detector shown has a plurality of 
sub-areas Ti to Tn, each of which comprises a plurality of image areas. In the present 
embodiment an image area is represented by one or more adjacently situated image columns 
which contain a plurality of adjacently situated sensor elements. The sub-area Ti contains 64 
5 image columns Si to S64. The adjoining sub-area T2 and all subsequent sub-areas T3 to Tn 
also contain 64 image columns each. With each sub-area Ti to Tn there is associated a 
respective read-out unit, which read-out units are in this case implemented as respective 
amplifiers Vi to Vn. 

The image data of the adjoining image columns 854 and Ses, amplified by the 
10 amplifiers Vj and V2, is applied to the analysis unit, that is, first to the histogram generator 15 
included therein. 

Fig. 3 shows the histograms of the two image columns 854 and Ses observed. It 
appears clearly that the image data of the image column 854 exhibits a maximum for a grey 
scale value around 13,500 and that the image data of the image column 855 has its maximum 

15 at y = 14,000. The histogram generator produces a histogram with a classification of the grey 
scale values in grey scale classes or areas, the frequency of occurrence of a grey scale value 
in the supplied image data of the image area observed being represented in the relevant class. 

Based on the assumption that the image data of adjoining image columns of 
neighboring amplifiers normally has the same grey scale values, the present case involves a 

20 difference which can be attributed to a difference in amplification behavior of the amplifiers 
Vi and V2. In the case of the same amplification behavior, the variations of the density 
functions f64(y) and fesiy) would register. 

Fig. 4 shows the integrated, or in this case the summed, variation of the 
density functions f64(y) and fesCy)- When the density functions f64(y) and fesCy) of the image 

25 columns 8^ and 855 are applied to the summing unit, the cumulative histograms or the 

distribution functions F64(y) and FesCy) of the density functions are formed. These cumulative 
histograms represent the yariation of the summing of the absolute grey scale values and 
indicate the probability that a grey scale value is below a given value. 

Fig. 5 shows the adaptation performed in the adaptation unit 17. To this end, 

30 the distribution function value F^siy) is determined from the distribution function (step 1) for 
each grey scale value yes. This distribution function value FgsCy) is used to find the 
distribution function value F64(y) (step 2). For this distribution function value F64(y) the grey 
scale value y64 of the distribution function F64(y) is determined. The adaptation table 
A(T2— >Ti) of grey scale values yes thus determined for grey scale values y64 represents a rule 
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in the form of a table for the adaptation of the image data of the sub-area T2 to the sub-area 
T,, 

Repeatedly performing these adaptation steps on each time two neighboring 
sub-areas yields correction data which is suitable for adapting these two sub-areas each time 



to the adaptation table A(T3^T2) in order to adapt the sub-area T3 to the sub-area T2 results 
in an adaptation table A(T3^Ti) for adapting the sub-area T3 to the sub-area Ti. Adaptation 
' tables are thus found for all amplification characteristics of the sub-areas Ti to Tn; these 
tables can be used directly in the correction unit with the same expenditure on time. 

Fig. 6 shows an X-ray image containing artefacts in the form of stripes. 

Fig. 7 shows the X-ray image corrected by means of the method according to 

the invention. 

A flat X-ray detector used has, for example, a surface with 1024 x 1024 sensor 
elements. Each grey scale value detected by this X-ray detector is encoded by way of 10 bits. 
The image repetition frequency amounts to approximately 30 images per second. Thus, it is 
necessary to refrain from applying each individual image to the analysis unit and to apply 
instead, for example, only every second image. The analysis unit calculates the correction 
data of the adaptation tables; for all single images to be processed during this calculation 
there is performed a correction by means of the respective instantaneous adaptation table 
which is stored, for example, in the memory. 



sub-area Ti to T5 in this detector again has 64 image columns. In order to make it 
unambiguously clear which pixel in the image belongs to which sub-area Ti to T5, that is, to 
optically visualize the boundaries G between the sub-areas Ti to T5, a given grey scale value 
is assigned to each pixel of a sub-area Ti to T5 in this image; this grey scale value is 
characteristic of the respective sub-area. Because the images of X-ray analysis devices or 
similar apparatus are customarily displayed in grey scale values, hereinafter image values are 
also taken as grey scale values. Evidently, however, the invention is not restricted to such 
grey scale values and any other image data may also be concerned; for example, in the case 
of color images it may be image data represented by a vector containing the RGB 
components. It may also concern images which are represented by amplitude and phase 
values. 

The method shown in Fig. 10a for the representation of the pixels, associated 
with a given amplification channel, by means of a given, predetermined grey scale value is 



in respect of their amplification characteristics. Application of the adaptation table A(T2— >Ti) 



Fig. 10a shows an image of a real detector with five sub-areas Ti to T5. Each 
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particularly useful to determine which pixel belongs to which amplification channel in the 
case of detectors having a complex structure. For the conventional detectors, being grouped 
in individual image columns and adjacently situated image columns being assigned a read- 
out unit, it suffices to define merely the horizontal co-ordinate constituting the boundary 
5 between two sub-areas. 

Fig. 10b shows an image of a hand which has been formed by means of the X- 
ray detector shown in Fig. 10a. The disturbing formation of stripes due to the different 
amplification behavior of the read-out units of the individual sub-areas Ti to T5 is clearly 
visible. 

10 Fig. 9 shows a device according to the invention which is intended to avoid 

such formation of stripes. The image produced by the detector 1 1 and prepared in a 
preprocessing unit 19 which performs, for example, an offset correction, is applied to a 
correction unit 13. The image values of the individual pixels of the individual sub-areas are 
corrected in the correction unit 13 in such a manner that the amplification curves of the read- 

15 out units Vi to Vn of the individual sub-areas Ti to Tn are adapted to one another. The 
corrected image is then displayed on a display unit 10. 

In order to acquire the correction values for the correction of the image data, 
the images from the preprocessing unit 19 are also applied to an analysis device 12. A 
respective correction factor KF is formed for the individual image values GW of a given sub- 

20 area Ti to Tn in the analysis device 12. These correction factors are then stored in the form of 
an adaptation table LUT in a memory 14 wherefrom they can be fetched at any time by the 
correction unit 13 in order to correct an image. Alternatively, of course, the analysis device 
12 can also apply the correction factors KF directly to the correction unit 13. 

In the present embodiment the correction factors KF are determined as shown 

25 in Fig. 8. In this Figure the grey scale values GW are plotted over the individual pixels along 
a row of the detector, that is, perpendicularly to a boundary G between two sub-areas Ti to 
Tn; in this case only the pixels Fei to Pe? of the last three image columns 852 to 854 of the first 
sub- area Ti and the first three image columns S^s to Se? of the second sub-area T2 are shown. 

As is shown in Fig. 8, the last pixel p64 of the first sub-area Ti has a given grey 

30 scale value GW64 and the first pixel Pas of the second sub-area T2, which pixel directly 

adjoins the pixel P64, has a different grey scale value GWes. This means that the grey scale 
value changes abruptly from the last pixel P64 of the first sub-area Tj, across the boundary G, 
to the first pixel Pes of the second sub-area T2. 
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For the correction method according to the invention it is assumed that this 
transient may be attributed to the different ampHfication characteristics of the two amplifiers 
Vi and V2. In order to match the different amplification characteristics, therefore, a search is 
made for correction values whereby the grey scale values of the second sub-area T2 are 
5 adapted to the grey scale values of the first sub-area Ti. 

To this end, first an estimated value SWgs is determined for the grey scale 
value GWes of the first pixel Pes of the second sub-area T2. This estimated value SWes is 
determined by extrapolation over the grey scale values GWes and GW64 of the two pixels Pes 
and P64 of the first sub-area Ti. The quotient of the estimated value SWes and the actual grey 
10 scale value GWes of the pixel Pes then yields a correction factor for all grey scale values in 
the sub-area T2 which correspond to grey scale value GWes of the pixel Pes observed. 
O In order to improve the adaptation, in the present embodiment an estimated 

SA value SWe4 is formed for the last pixel P64 of the first sub-area Ti in the same way by an 

[5 analogous extrapolation at the boundary G from the right-hand side, that is, over the grey 

ffi 15 scale values GWes and GWee of the pixels Pes and Pee- A second correction factor is derived 
LI from the ratio of the actual grey scale value GW64 to the estimated value SWe4; this second 

^ correction factor corresponds to the first correction factor in the ideal case. A mean value can 

fy be formed from these two correction factors; this mean value is ultimately used as the 

correction value determined for the pixel Pes- 
S 20 This method is performed for all pixels of the first image column Ses of the 

second sub-area T2. Generally speaking, a plurality of correction factors will then be 
determined for a given grey scale value GW, because this grey scale value GW may possibly 
be present at various pixels of the image column Ses- These correction values for each time 
the same grey scale values are accumulated and, for example, a mean value is formed 
25 therefrom. Evidently, it is also possible, and possibly even better, be it also more calculation 
time consuming, to use other means for combining the various correction values for one and 
the same grey scale value of a sub-area. For example, all correction values for a grey scale 
value can be taken up in a histogram, after which the modulus or the median is determined. 
Individual statistical errors of the method are thus corrected particularly well. 
30 This procedure within the analysis device 12 is shown again in Fig. 9. In a first 

unit 20 two estimated values SWes and SW64 are formed from the right and from the left by 
extrapolation at a boundary G. The unit 21 then forms a correction factor KF(GW) for a 
given grey scale value GW from these estimated values SWes and SW64. In a subsequent unit 
22 all correction factors KF are averaged for a given grey scale value GW so as to obtain the 
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ultimate correction factor for this grey scale value GW. The analysis device 12 of the 
memory 14 and the correction unit 13 may in principle be formed by a conMnon computer; 
the individual units 20, 21, 22 could then also be implemented as appropriate software in the 
computer. 

This procedure for determining the correction factors is performed for each 
sub-area, a sub-area Tj being adapted each time to the directly preceding, neighboring sub- 
area Tj-i. This means that the grey scale values of the sub-area T2 are adapted to the sub-area 
Ti, the grey scale values of the sub-area T3 to the sub-area T2 etc. All sub-areas T2 to Tn are 
thus successively adapted to the sub-area Ti. 

In order to obtain a respective absolute correction value, related to the first 
sub-area Ti, for each sub-area T2 to Tn, the relative correction values of the individual sub- 
areas Tj (related to their preceding neighboring sub-area Ti.i) must be multiplied each time by 
the absolute correction values Ti.i of the preceding neighboring sub-area. 

It must be emphasized again that this algorithm does not eliminate the non- 
linearity of the individual amplifiers, but only adapts the non-linearities of the channels to 
one another. The correction factors thus determined for the individual grey scale values for 
each individual sub-area Ti to Tn are then taken up in an adaptation table LUT in the memory 
14. 

Normally speaking, not all feasible grey scale values GW necessarily occur in 
the respective first image columns of a sub-area Ti to Tn, because most detectors involve 
several tens of thousands of different grey scale values. The adaptation table can then simply 
remain empty for these grey scale values and interpolation methods can be used so as to fill 
in the values. The size of the adaptation table can also be reduced, for example by combining 
the correction values for a plurality of grey scale values in the adaptation table. That is, a 
plurality of grey scale values are combined and provided with a common correction factor. 
This makes sense notably when the differences in non-linearity vary only slowly as a 
function of the grey scale value. 

Fig. 10c shows an adaptation table LUT for the five sub-areas Ti to T5 of the 
detector shown in the Figs. 10a and 10b. A grey scale value is plotted each time as a 
correction factor which replaces the respective grey scale value occurring in the sub-area. In 
the representation shown the grey scale values GW to be replaced for the sub-areas Ti to T5 
increase from the top downwards in the marked direction of the arrow. For each sub- area Ti 
to T5 the adaptation table LUT contains only a single vector of entries which extends in the 
direction of the arrow. The expansion in the width is performed merely for a better 
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representation of the different sub-areas Ti to T5. As is clearly shown, the grey scale values 
themselves are entered only in the sub-area Ti. This means that the sub-area Ti is corrected 
quasi with a correction factor 1.0, because all further sub-areas T2 to T5 are adapted to this 
first sub-area Ti. 

Finally, Fig. lOd shows the X-ray image of Fig. 10b after application of the 
correction values from the adaptation table LUT of Fig. 10c to the image data. It is clearly 
shown that the undesirable stripes have been fully eliminated. 

An extrapolation over each time two pixels was performed so as to find the 
estimated value in the illustrated example. It will be evident that an extrapolation over several 
adjacently situated pixels can also be performed. It is, of course, also possible to perform the 
method without such extrapolation. For example, in the simplest case direct adaptation can be 
performed between two neighboring pixels. This means that for the estimated value for the 
grey scale value of a given pixel in the first image column 855 of the second sub-area T2 use 
is made of the grey scale value of the neighboring pixel in the last image column 854 of the 
first sub-area Ti, etc. 

Moreover, of course, this method is not restricted to the use for detectors with 
only vertically subdivided sub-areas Ti to Tn- Evidently it can also be used for detectors 
which are subdivided, for example into a plurality of vertical sub-areas and additionally in a 
plurality of horizontal sub-areas. 

Various methods can be applied in order to enhance the reliability of the 
method and to avoid errors due to natural gradients in the image which are not related to the 
boundary G between two sub-areas Ti to Tn; they can also be used in combined form: 

1. Only correction values of the pixels along a boundary for which the correction 
values do not deviate excessively from one another because of the right-hand and left-hand 
extrapolation are taken into account. 

2. Correction factors which lie outside a given, predetermined correction factor 
range are not taken into account. 

3. Correction values for which the gradient of the grey scale value relative to the 
grey scale value (that is, the slope of the gradient) exceeds a maximum value fixed between 
the sub-areas (absolute limit value of the gradient) are not taken into account. 

4. Correction values which deviate too much from the correction values of a 
previously determined adaptation table are not taken into account, 

5. A smoothing operation is performed on the non-empty entries in the 
adaptation table so as to allow only small variations of the correction factors. 
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6. All entries of the adaptation table are reduced as soon as a given maximum 

deviation value is exceeded. This restriction is related to the maximum non-linearity 
specification of the individual amplification channels and can best be illustrated on the basis 
of an example. For example, when the correction factor within an adaptation table of a given 
5 sub-area amounts to 1.01 in the case of a low grey scale value and to 1.05 in the case of a 
high grey scale value, the non-linearity of this sub-area, in comparison with the neighboring 
sub-area, amounts to (1.05 - 1.01) x 100% = 4%. However, when it is known on the basis of 
the specification of the amplifier that the non-linearity amounts to only 0.5%, the maximum 
non-linearity between two channels can inherently amount to no more than twice the non- 
10 linearity of a single amplifier only, that is, 1%. In this case the non-linearity entered in the 
adaptation table can be reduced from 4% to 1% by dividing the entire range by 4. More 
y specifically, in as far as the entries in the adaptation table concern correction factors, the 

SI value 1 can then first be subtracted from each individual entry, after which the overall mean 

^ value is determined for all entries and the deviations from the mean value are divided by the 

J 15 factor 4. It will be clear that such an adaptation requires an extremely large amount of 
M: calculation work. 

The method can be used for single exposures in that the correction values 
nJ determined from the image are applied to the image itself. However, it can be used equally 

11 well for a sequence of successive images, in which case the correction values are 

y 20 advantageously derived from an adaptation table which is formed from the image data of the 
preceding images and is checked and corrected at regular intervals by way of correction 
values derived from newly acquired images. 

The method can be used for a wide range of different image pick-up systems. 
It serves mainly to eliminate the visibility of the differences in the non-linearity of the 
25 preamplifiers in flat dynamic X-ray detectors (FDXD), but is not restricted to such an 

application. The XTV16 image sensor is a further example of an image pick-up system in 
which the method can be used. In this detector different sensors and preamplifiers are used 
for the even and the odd image columns. Because each sub-area in this detector consists of 
only a single column, however, the method of extrapolation over several pixels cannot be 
30 used; adaptation is then performed by comparison of neighboring pixels. Most of the methods 
described for the above embodiment for the combination of different correction values of one 
and the same grey scale value and for the enhancement of the reliability can also be 
performed for such a comparison of individual pixels without extrapolation. 
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The method can in principle be used even for performing a masking correction 
by defining a set of concentric circles as different sub-areas of the individually arriving 
image. 



